* Stata statistical program accompanying...
* Edward D. Mansfield and Eric Reinhardt
* "Multilateral Determinants of Regionalism:
* The Effects of GATT/WTO on the Formation of Preferential Trading Arrangements"
* International Organization 57 (Fall 2003), pp. 829-862

* For questions, contact:
* Eric Reinhardt
* Department of Political Science
* Emory University
* Atlanta, GA 30322
* erein@emory.edu
* 404-727-4977
* http://userwww.service.emory.edu/~erein/

* This syntax and output is for/by Stata SE version 8.2.

set more off
set matsize 1000
set mem 325m
use pta2, clear
log using mansfield-reinhardt-2003.log, append

**************************************************************
**************************************************************
* Published Models *
**************************************************************
**************************************************************

******* MODEL 1 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)

******* MODEL 4 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)

**************************************************************
**************************************************************
* Various statistics reported in the paper *
**************************************************************
**************************************************************

* notes 89 and 91, page 848 *
******* MODEL 1 *******
quietly logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
testparm gattl1 round2 bad14 f_14
testparm spldy*
******* MODEL 2 *******
quietly logit newptar distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
testparm gattl1 round2 bad14 f_14
testparm spldy*
******* MODEL 3 *******
quietly logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)
testparm dtgattl1 round2 bad14 f_14
testparm spldy*
******* MODEL 4 *******
quietly logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
testparm gattl1 round2 bad14 f_14
testparm spldy*
******* MODEL 5 *******
quietly logit newptar distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
testparm gattl1 round2 bad14 f_14
testparm spldy*
******* MODEL 6 *******
quietly logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
testparm dtgattl1 round2 bad14 f_14
testparm spldy*

* note 96, page 852: what if set missing trade values to zero? *
replace tr95=0 if tr95==.
******* MODEL 1 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 4 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
use pta2, clear

* tests in bottom paragraph page 851, top page 852 *
* Country Fixed-Effects *
tab stnum, gen(fes1_)
tab stnum2, gen(fes2_)

******* MODEL 1 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* fes1_* fes2_* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* fes1_* fes2_* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* fes1_* fes2_* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 4 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* fes1_* fes2_* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* fes1_* fes2_* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* fes1_* fes2_* if pta2rl1~=1 & side_a==n, robust cluster(dyad)

use pta2, clear
* tests in bottom paragraph page 851, top page 852 *
* control for hegemony, MID, trade dependence *
* 1st version of hegemony variable *
******* MODEL 1 *******
logit newptar hegl1 tradep1 mid15 alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar hegl1 tradep1 mid15 distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar hegl1 tradep1 mid15 alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 4 *******
logit newptar hegl1 tradep1 tradep2 mid15 alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar hegl1 tradep1 tradep2 mid15 distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar hegl1 tradep1 tradep2 mid15 alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
* 2nd version of hegemony variable *
******* MODEL 1 *******
logit newptar top2rat tradep1 mid15 alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar top2rat tradep1 mid15 distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar top2rat tradep1 mid15 alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 4 *******
logit newptar top2rat tradep1 tradep2 mid15 alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar top2rat tradep1 tradep2 mid15 distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar top2rat tradep1 tradep2 mid15 alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)

* tests in bottom paragraph page 851, top page 852 *
* drop dyads with US in them *
use pta2, clear
drop if stnum==2 | stnum2==2
******* MODEL 1 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 4 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)

* tests in bottom paragraph page 851, top page 852 *
* drop dyads with at least one EC founding member *
use pta2, clear
gen ec61=1 if country1=="Belgium" | country1=="Netherlands" | country1=="France" | country1=="Luxembourg" | country1=="Germany" | country1=="Italy"
gen ec62=1 if country2=="Belgium" | country2=="Netherlands" | country2=="France" | country2=="Luxembourg" | country2=="Germany" | country2=="Italy"
drop if ec61==1 | ec62==1
******* MODEL 1 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 2 *******
logit newptar distance demp981 g1 gc1 pdsrfc pdsrfc2 tp2r1 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 3 *******
logit newptar alliance distance demp981 tr95 g1 gc1 growth1 pdsrfc pdsrfc2 tp2r1 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1, robust cluster(dyad)
******* MODEL 4 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 5 *******
logit newptar distance demp981 demp982 g1 g2 gc1 gc2 pdsrfc pdsrfc2 tp2r1 tp2r2 year gattl1 round2 bad14 f123 f_14 spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)
******* MODEL 6 *******
logit newptar alliance distance demp981 demp982 tr95 g1 g2 gc1 gc2 growth1 growth2 pdsrfc pdsrfc2 tp2r1 tp2r2 year dtgattl1 round2 bad14 f123 f_14 fmrcolb postcom spldy* if pta2rl1~=1 & side_a==n, robust cluster(dyad)

log close

exit, STATA clear
